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ABSTRACT 



Torque teno viruses (TTVs) are recently dis- 
covered DNA viruses, with heterogeneous genomes, 
highly prevalent in populations worldwide. The spe- 
cies that infect humans are Torque teno virus (TTV), 
Torque teno midi virus (TTMDV) and Torque teno 
mini virus (TTMV). High-resolution melting analysis 
(HRMA) is a sensitive and effective method for ge- 
notyping and mutation scanning. Up to now, HRMA 
has not been utilized for detection of TTVs. 

The aim of this study was to asses if HRMA 
is suitable for detecting TTVs variants. DNA was 
extracted from the blood and saliva of 13 healthy 
subjects for method optimization. Additionally, saliva 
samples from 1 00 healthy individuals were collected 
for estimating the TTVs' prevalence. Viral DNA was 
amplified by heminested polymerase chain reaction 
(PCR). Second round amplicons were used for the 
HRMA. The samples were analyzed using two fluo- 
rescent dyes, SYBR® Green I and EvaGreen®. 

The prevalence values for TTV, TTMDV and 
TTMV were 7 1 .0, 3 1 .0 and 54.0%, respectively. The 
three major melting curve patterns corresponding to 
TTV, TTMDV and TTMV on HRMA can be eas- 
ily distinguished regardless of kit used. Our results 
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showed that HRMA is a rapid and efficient method 
of detecting human TTVs. 

Keywords: Torque teno virus (TTV); Torque 
teno midi virus (TTMDV); Torque teno mini virus 
(TTMV); High-resolution melting analysis (HRMA). 

INTRODUCTION 



Torque teno viruses (TTVs) are small, non envel- 
oped viruses, with an icosahedral capsid and a single- 
stranded circular DNA genome. One characteristic 
of TTVs is the extreme genomic heterogeneity that 
is highly uncommon for DNA viruses. Three virus 
species that infect humans have now been identified: 
TTV {Torque teno virus) [1], TTMV {Torque teno 
mini virus) [2] and TTMDV {Torque teno midi virus) 
[3]. Torque teno mini virus (2.8-2.9 kb) and TTMDV 
(3.2 kb) have smaller genomes than the TTV genome 
(3.75-3.9 kb). The TTV, TTMV and TTMDV each 
consist of a wide repertoire of isolates with genome- 
wide sequence divergences of at least 50.0, 40.0 and 
33.0%, respectively [4]. Based on genomic hetero- 
geneity, TTVs were separated into 39 genotypes and 
five major genetic groups [5,6], and it was suggested 
that these viruses may exist as a swarm of closely 
related but different viral species [7]. 

The sequence divergence is not evenly distrib- 
uted throughout the genome. The untranslated region 
(UTR) is well conserved (73.0% nucleotide identity 
between the two highly divergent isolates TA278 
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Table 1. Amplification conditions and primers [4] for heminested polymerase chain reaction. 



First Round of Polymerase Chain Reaction 




90°C 


55°C 


72°C 


1 x a 


10 min. 






1 X 


2 min. 






2 x 


2 min. 


30 sec. 


30 sec. 


35 x 


35 sec. 


30 sec. 


30 sec. 


1 X 






1 min. 


Second Round of Polymerase Chain Reaction 


1 x a 


10 min. 






1 X 


2 min. 


30 sec. 


30 sec. 


27 x 


30 sec. 


30 sec. 


30 sec. 


1 X 






1 min. 


Primers 


Primer Code 


PCR Round 


Specificity 


Forward of Reverse 


NG779 


first and second (TTV) 


general 


F 


NG780 


first and second (TTV) 


general 


F 


NG781 


first 


general 


R 


NG782 


first 


general 


R 


NG785 


second 


TTV 


R 


NG795 


second 


TTMDV 


F 


NG796 


second 


TTMDV 


R 


NG791 


second 


TTMV 


R 


NG792 


second 


TTMV 


F 


NG793 


second 


TTMV 


F 


NG794 


second 


TTMV 


F 



n Additional step for hot-start enzymes (Immolate™ DNA polymerase, Bioline; and Dry PCR Master Mix, Bioron). 



and SANBAN vs. 57.0% for the entire genome) and 
contains several highly con-served sequences (with 
more than 90.0% identity between isolates). In con- 
trast, the translated region is characterized by a very 
high degree of diversity [8,9]. 

The TTVs disease-causing potential is still under 
debate. Rocchi et al. [ 1 0] showed that certain variants 
pre-sent CpG motifs that are activators of proinflam- 
matory cytokine production via toll-like receptor 9, 
and thus, may increase severity of inflammatory dis- 
eases. For the identification of these variants, new 
methods for assessing the viral diversity are needed. 

Polymerase chain reaction (PCR)-based methods 
were used to detect the presence of TTVs DNA in bio- 
logical samples. The high TTVs sequence divergence 



and the particularities of PCR assay used for detec- 
tion (primers, genomic region amplified, protocols) 
[11] may determine preferential amplification of some 
TTV genotypes and, thus, can influence the preva- 
lence rates of TTVs presented in different articles. 

One of the first diagnostic assays for TTVs was 
based on PCR that used primers derived from the N22 
region of ORF1 [1]. This primer set did not amplify 
all TTV phylo-genetic groups. Subsequently, PCR 
techniques aimed for the conserved UTR, near the 
5' end of ORF2 (Takahashi et al. [12]) or near the 3' 
end of ORF3 (Leary et al. [13]). 

Methods based on analysis of changes in fluo- 
rescence intensity contribute to genome diversity 
analysis. One of these methods is high-resolution 
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Figure 1. Different patterns of melting curves (a) and the difference graphs (b) for TTV (red), TTMDV (blue) 
and TTMV (green) in HRMA carried out with SensiMix ® EvaGreen® fluorescent dye (Bioline). 



melting analysis (HRMA) [14], which uses high data- 
density acquisition during melting of amplicons in the 
presence of intercalating dye. Because of its low-cost, 
speed and simplicity, this method is used for genetic 
discrimination of various microorganisms, includ- 
ing viruses [15-18] and bacteria [19,20]. However, 
as far as we know, it has never been used for TTVs 
genotyping. The aim of this study was to assess if 
HRMA can be used for efficient and rapid detection 
of TTVs variants. 

MATERIALS AND METHODS 

Samples and DNA Extraction. A total of 113 
healthy subjects from Bucharest, Romania, ranging in 
age from 19 to29 years, were selected for this study. 
Blood and saliva samples from 13 subjects were used 
for protocol optimization, whereas saliva samples 
(n = 100) were used to estimate the TTVs preva- 
lence. All subjects gave written informed consent 
to participate in this study according to the Helsinki 



Declaration of 1975, as revised in 1983, concerning 
the ethical principles for medical research involving 
human subjects. 

DNA was extracted from each blood sample 
using four different commercial kits (recommend- 
ed for DNA extraction from eukaryotic cells): two 
spin-column-based protocols, AxyPrep™ Blood 
Genomic DNA Miniprep Kit (Axygen, Union City, 
CA, USA) and Ron's Blood DNA Mini Kit (Bioron, 
Ludwigshafen, Germany), and two pro-tocols without 
spin-column. Wizard® Genomic DNA Purification 
Kit (Promega, Madison, WI, USA), Pronto® DNA 
Extraction Kit (Pronto Diagnostics, Rehovot, Israel). 
The DNA was extracted according to the manufac- 
turer's instructions from 250 uL of blood. A one-tube 
purification method, OTP reagent (Bioron), was used 
for DNA extrac-tion from saliva. 

Polymerase Chain Reaction. The detection of 
the three TTVs was realized using a two-step hemi- 
nested-PCR protocol as described by Ninomiya et 
al. [4]. Each DNA sample was PCR amplified using 
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Table 2. Torque teno virus genotyping results obtained with three DNA polymerases from blood and saliva DNA 
extracted with different commercial kits. (A: GoTaq® DNA Polymerase, Promega; B: Immolate™ DNA 
Polymerase, Bioline; C: Dry PCR Master Mix, Bioron) 





DNA Extraction Kit 




AxyPrep™ Blood 
Genomic DNA 
Miniprep Kit 


Ron's Blood 
DNA Mini Kit 


Wizard® Gnomic 
DNA Purification 
Kit 


Pronto® DNA 
Extraction Kit 


OTP Reagent 
(saliva) 




DNA Polymerase 


n 


A 


B 


C 


A 


B 


C 


A 


B 


C 


A 


B 


C 


A 


B 


C 


1 






















[-] 










1 


[+] 


[+] 


[+] a 


[+] a 


[+] 


[+] a 


[+] 


[+] 


[+] 


[+] a 


[+]• 


[+] 


[+] a 


[+] a 


[+] a 


5 


[+] 


[+] a 


[+] a 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


4 


[+] a 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


5 


[+] 


[+] a 


[+] 


[+] 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+]' 


[+] 


[+] 


[+] 


[+] 


6 


[+] 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] a 


7 


[+] 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


8 






















[-] 










9 


[+] a 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+]' 
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[+] 
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[+] a 
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12 
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[+] a 
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[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] 


[+] a 
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13 
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[+] 


[+] 
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Indicates a weaker signal in agarose ge. 




Figure 2. Melting curve patterns for HRMA carried out Figure 3. Melting curve patterns for HRMA carried out 



using SYBR Green Master Mix (Invitek). using Maxima® SYBR Green qPCR (Fermentas). 
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Table 3. Amplicon properties for five isolates of the tor- 
que teno virus, torque teno midi virus and torque teno 
mini virus. 



Isolate 


Length 
(bp) 


%GC 


T 

m 


Torque Teno Virus 


TA278 


114 


63 


85 


JA9 


114 


63 


85 


JA1 


114 


65 


86 


JA4 


114 


65 


86 


JA20 


114 


66 


86 


Torque Teno Midi Virus 


MDJHem2 


88 


67 


85 


SAV2 


88 


68 


85 


MD2-032 


88 


69 


86 


MD 1-073 


88 


70 


86 


SAV1 


88 


70 


86 


Torque Teno Mini Virus 


CBD203 


70 


54 


77 


CLC062 


72 


58 


79 


CBD279 


72 


60 


80 


CBD231 


72 


61 


80 


CLC138 


72 


61 


81 



three DNA-polymerases: GoTaq® DNA Polymerase 
(Promega), Immolase™ DNA Polymerase (Bioline, 
Luckenwalde, Germany), and Dry PCR Master Mix 
(blue-stain) (Bioron). 

For the first step, the amplification reaction was 
car-ried out in a total volume of 1 5 uL containing 1 x 
PCR buffer solution, 1.5 mmol/LMgCl 2 , 100 umol/L 
dNTP mix, 0.5 umol/L of each primer, 0.6 units of 
DNA-polymerase and 1.5 uL of genomic DNA ex- 



tract. The second step consisted of three different 
reactions for the specific amplification of the viral spe- 
cies. For the second step, 1 uL of amplicon from the 
first round of PCR was used as a template. The ampli- 
fication programs and primers can be found in Table 
1 . The amplicons of the second PCR were resolved 
by 2.0% (w/v) agarose gel and were electropho-resed 
at 5 V/cm for 20 min. The gel was visualized under 
UV light after ethidium bromide staining. 

High-Resolution Melting Analysis. The sam- 
ples co-infected with all three TTVs were further 
analyzed by HRMA using three commercial kits 
and two different fluorescent dyes: SensiMix® with 
EvaGreen® (Bioline), Invitek SYBR Green Master 
Mix (Invitek, Berlin, Germany) and Maxima® SYBR 
Green qPCR (Fermentas, Glen Burnie, MD, USA). 
The amplification reaction was carried out in a total 
volume of 20 uL and used 2 uL of second-round 
amplicons as template DNA. High-resolution melt- 
ing analysis was performed in a Rotor-Gene 6000 
series instrument (Corbett Research, Sydney, NSW, 
Australia) using software version 1.7. 

Melting curves were generated by increasing 
the tem-perature from 65°C to 90°C ramping by 0. 1 
degrees per step, with a gap of 2 seconds between 
steps. For the HRMA normalization regions of 75- 
77°C and 88-90°C were applied. 

RESULTS 



The prevalence of TTV, TTMDV and TTMV in 
saliva samples was 7 1 .0, 3 1 .0 and 54.0%, respectively. 
There were no qualitative differences in the genotyping 
results obtained for the same sample extracted with 
four commercial kits and amplified with three different 
DNA polymerases. The distribution of TTVs did not 
differ in blood and saliva samples (Table 2). 




Figure 4. The HRM amplicons electrophoresed in 2.0% (w/v) agarose gel stained with 
ethidiumbromide.Lanes2,5and8:TTV;lanes3,6and9:TTMDV;lanes4,7andl0:TTMV; 
lane 1: pUC19 DNA/MspI (Hpall) marker (Fermentas). 
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Ten samples with triple coinfection were selected 
for the HRMA using SensiMix®, SYBR Green Master 
Mix and Maxima® SYBR Green qPCR commercial kits 
(Figure 1 a, Figure 2 and Figure 3). Prior to HRMA, am- 
plicons were verified by gel electrophoresis (Figure 4). 

The three major melting curve patterns corre- 
sponding to TTV, TTMDV and TTMV can be easily 
distinguished regardless of the kit used, although they 
present slightly different aspects. The melting curve 
patterns were similar when samples were retested 
under the same conditions. 

DISCUSSION 



Torque teno viruses have several characteris- 
tics (e.g., genomic heterogeneity, great number of 
isolates, high prevalence) that play an important 
role in the way their analysis must be approached. A 
wide range of TTV prevalence has been described 
worldwide. Depending on the identification method 
used, TTV prevalence may vary from 46.0-62.0% in 
Brazil [21,22] to 94.0% in Russia [23]. For TTMV, 
prevalence ranges between 48.0% in Norway [24] 
and 67.0-72.0% in Brazil [21,22], and TTMDV oc- 
curs in at least 40.0% of the general population [25]. 
Our results are in accordance with prevalence values 
reported for TTVs in other populations, except for 
TTMDV, which has a lower prevalence in subjects 
selected for this study (31.0%). 

Analyzing such heterogeneous genomes rep- 
resents a challenge. The high sequence variation of 
the genomes, along with high percentage of cyto- 
sine and guanine, make primer design difficult. The 
heminested PCR assay designed by Ninomiya et al. 
[4] amplifies at least 49 TTV isolates, 20 TTMDV 
isolates, and 13 TTMV isolates. 

The 1 0 melting curve patterns on the HRMA sug- 
gest that different viral isolates may be present in our 
cohort, or individuals may be coinfected with several 
isolates. These may explain the shift of the melting 
curve for TTMDV and TTMV samples marked in 
Figure 1 and Figure 3. 

The melting curve aspect is influenced by the 
length, percentage of cytosine and guanine (%GC) 
and melting point (T) of the amplification products. 
The %GC and T m were calculated in silico for ampli- 
cons produced by five isolates of each virus (Table 2). 

The differences between TTV and TTMDV am- 
plicons are minimal, and an accurate discrimination 



is difficult. Moreover, we have run a folding simula- 
tion for each of the amplicons using Quickfold [26] at 
65°C (the temperature in the beginning of HRMA) and 
observed that each amplicon forms slightly different 
secondary structures when denatured. This influences 
the melting curve pattern. Subsequent to HRMA, the 
resulting amplicons were verified by gel electropho- 
resis and no differences were observed (Figure 4). 

Another component of the HRMA that influ- 
ences the curve patterns is the fluorescent dye. There 
are various types of double-stranded DNA (dsDNA) 
intercalating dyes with different properties. For 
HRMA, the dye must provide detailed information on 
the melting behavior of an amplified target. Ideally, 
the dye should not bind preferentially to pyrimidines 
or purines, change the T m of the amplicon, or inhibit 
DNA amplification. 

There are two main types of dyes: saturating 
and non saturating dyes. SYBR® Green I is a non 
saturating dsDNA intercalating dye and is not usu- 
ally recommended for high-resolution melt applica- 
tions because at high concentrations, SYBR® Green 
I inhibits the DNA polymerase. At low concentra- 
tions SYBR® Green I is able to redistribute from the 
melted regions back to the regions of dsDNA, which 
results in poor base-difference discrimination [27]. 

EvaGreen® is a special kind of saturating dye, 
so - called "release-on-demand." For this dye, the fluo- 
rescence is quenched when unbound to DNA, this al- 
lows the use of non saturating dye concentrations, thus 
ensuring no PCR inhibition. EvaGreen® improves 
the resolution and accuracy of HRMA by increased 
fluorescence and lack of redistribution during melting 
[28]. Farrar et al. [29] showed that EvaGreen® is more 
suitable than SYBR® Green I for HRMA. 

Our results show that the best discrimination is 
obtained using EvaGreen® (Figure 1 vs. Figures 2 
and 3). The differences between the curve patterns 
obtained with SYBR® Green I may be due to the 
DNA-polymerase activity [SYBR Green Master Mix 
(Figure 2) vs. Maxima® SYBR Green qPCR (Figure 
3)] and/or dye concentration in the master mixes used 
for this study. 

As with any technique, HRMA analysis has its 
limitations. Fluorescent dyes used in HRMA lack 
sequence specificity and can bind to any dsDNA, 
including non targets such as primer dimmers and 
non specific products, which will bias the results of 
melting analysis [28]. In addition, all the PCR com- 
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ponents are present at the time of melting analysis 
and may have a great impact on the melting curve 
shape and position [30]. 

The amplicons corresponding to TTV, TTMDV, 
and TTMV obtained in the assay we used have dif- 
ferent lengths (Table 3). The %GC and also GC dis- 
tribution in relation to the ends/center of amplicons 
differ due to genomic heterogeneity These aspects 
influence the melting curves' aspect as well. 

Despite these limitations, HRMA is a sensitive 
method and provides more information on the amplifi- 
cation products, such as sequence-dependent shape of 
the melting curve and T , enabling discrimination of 
products with same length but different sequence [31]. 

Our results showed that HRMA is a rapid method 
of detecting human TTVs (HRMA takes approxi- 
mately 20 min. after second round amplicons are ob- 
tained) compared to the classical PCR-electrophoresis 
method, which is more time-consuming (gel prepara- 
tion, running and staining). High-resolution melting 
analysis provides additional information regarding 
amplification products (T m , melting curve shape) 
compared to classic PCR methods followed by gel 
electrophoresis, which indicate only the presence or 
absence of the target sequence. 

In conclusion, due to the advantages of this tech- 
nique, HRMA is a rapid and accurate method for 
detecting TTVs. Developing new and more sensitive 
HRMA assays may lead to easy and accurate detec- 
tion of TTV isolates. 
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